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Abstract 

We introduce the neyative binomial states with negative binomial distribution as 
their photon number distribution. They reduce to the ordinary coherent states and 
Susskind-Glogower phase states in different limits. The ladder and displacement op- 
erator formalisms are found and they are essentially the Perelomov's su(l, 1) coherent 
states via its Holstein-Primakoff realisation. These states exhibit strong squeezing ef- 
fect and they obey the super-Poissonian statistics. A method to generate these states 
is proposed. 

PACS: 03.65.-w, 05.30.ch, 42.50.Ar 



*On leave of absence from Institute of Theoretical Physics, Northeast Normal University, Changchun 
130024, P.R.China. E-mail: hcfu@yukawa.kyoto-u.ac.jp 



1 Introduction 



Nonclassical states of the radiation fields, such as the number states, the coherent states and 
the phase states, play important roles in quantum optics and are extensively studied |1|]. The 
binomial states (BS) introduced by Stoler et al. in 1985 f2|, interpolate between the most 
nonclassical number states and the most classical coherent states, and reduce to them in two 
different limits. Some of their properties P, [| |J, methods of generation p], |], H, as well as 
their interaction with atoms 0, have been investigated in the literature. The notion of BS 
was also generalised to the intermediate number- squeezed states |7|, |J and the number-phase 
states f9j, the hypergeometric states as well as their g-deformation |JTT . 



The photon number distribution of BS is the binomial distribution in the probability 
theory ]12|, |I3fl . In this letter, we shall introduce and study negative binomial states (NBS) 
whose photon distribution is the negative binomial distribution (NBD) Different 



from the BS, the NBS are the intermediate phase- coherent states in the sense that they 
reduce to the Susskind-Glogower (SG) phase states |TJ]] and coherent states in two different 



limits (Sec. 2). We also derive their ladder and displacement operator formalisms and find 
that they are essentially the Perelomov's su(l, 1) coherent states via its Holstein-Primakoff 
(HP) realizations (Sec. 3). The NBS exhibit strong squeezing effects, but are not of sub- 
Poissonian statistics (Sec. 4). A method to generate these NBS is proposed in Sec. 5 and a 
summary with special emphasis on the comparison with BS is given in Sec. 6. 

2 Negative binomial states and their limits 

We define the NBS as 



ne* 9 ; M)~ = £ ^/fl-fo; M)e me \\n), (2.1) 



n=0 



in which {||n) \n — 0, 1, . . .} are the number states of the single-mode radiation field 

[M f ] = l, 6||0> = 0, ||n> = ^||0>, (2.2) 



M is a fixed positive integer, rf is the probability satisfying < rf < 1, and 

B-( m M)=^ M+ ^- i y^l-r ] Y I , n = 0,l,.... (2.3) 
The B~(rj; M) is called the NBD [F3] since it can also be written as 

B-( V - 1 M)^(l-rj 2 ) M (~^)(-V 2 ) n , n = 0,l,..., (2.4) 
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which has the similar form as the binomial distribution except for the two minuses and that 
n runs to infinity. The states ( |2.1| ) are referred to as the NBS since their photon distribution 
\(n\\rie t9 ;M)~\ 2 = B'faM) is the NBD. 

The parameter 9 (0 < 9 < 2tt) has clear physical meaning: it reflects the time de- 
velopment of the NBS. This can be seen from e~ lHt \r]e ie ; M)~ = Ir/e^ 6 '^; M)~ , where 
H = uj(N + 1/2) is the Hamiltonian of the single mode radiation field. 



As a probability distribution B n (77; M) satisfies [12 



XXfa;Af) = l, (2-5) 

n=0 

which means that the NBS are normalized. 
Let us consider two limiting cases: 

(1) . In the limit 77 — ^ 0, B~(rj; M) — > 5 n0 and thus the NBS go to the vacuum state. 

(2) . When 77 -> 0, M -> 00 with fixed finite r? 2 M = a 2 , the NBD goes to the Poisson 



distribution B~(r);M) —> e~ a a • jn\ |12]| . Accordingly, the NBS degenerate to the ordinary 
coherent states. 

We know that the BS degenerate to the number state in a certain limit ||. However the 
NBS do not maintain this feature. Instead, the NBS tend to the SG phase states in a certain 
limit. To achieve this, let us consider the case M — 1. In this case, \r\e %6 \ M)~ is simplified 

as 



V e ie ;l)- = Jl -r] 2 J2v n e m6 \\n). (2.6) 



n=0 



The photon distribution | (n\\i]e' 18 ; 1)~| 2 in this case is (1 — t] 2 )r] 2n , the geometric distribution. 
For this reason we call \r]e ie ; 1)~ the geometric states. One can easily verify that the geometric 
state \r]e ld ; 1)~ is the eigenstate of the SG phase operator E~ jnj with the eigenvalue r\e ld 

tf- = X»(n+l||, (2.7) 

n=0 



which is related to the annihilation operator b through polar decomposition b = E y N 



Now, multiplying a constant l/y27r(l — rf) to 1 77c ; 1)~ and then taking the limit 77 — > 1, 
we obtain the phase states 

1 1 00 

\9) = lim \ V e w ; 1)" = _ £ e ^||n), IT |0) = e J ». (2.8) 

^ x */27r(l - r/ 2 ) V2vr ^ 



The |#) states are non-normalisable, nonorthogonal, but resolve the identity 

d0|0)(0| = l. (2.9) 



The phase-state representation based on ( [2.9|) is a useful calculational tool [HJ . 

In this sense, we find that the NBS are the intermediate phase- coherent states. This is 
an important feature of the NBS. 



3 Displacement and ladder operator formalisms 



Let us recapitulate the displacement operator formalism of NBS |L3[ using the identity method 
developed in a previous paper |]16| . To this end, let us rewrite NBS ( j2.1|) (r/c = r/e zd ) 



foe; M)" = (1 - to c | 2 )^ £ — — -(^H^lO). (3.1) 



n=0 ™ ! 



Then, by making use of the following identity 

(tfg(N)) n \\Q) = (tf) n g(Q)g(l)---g(n-l)\\0), with g(N) = VWTN, N = 6+6. (3.2) 
we can write Eq. (|3.1| ) in the exponential form 



M 



\r, c - M)~ = (1 — tocl 2 )^ exp [r?c K+] ||0), 
/C+ = bWM + N = VM - 1 + iV6+. (3.3) 

It is interesting that /C + along with 

/C_ = (/C+) f = VM + iV6 = 6VM + AT — 1, /C = y + iV (3.4) 

generates the su(l, 1) algebra via its HP realization with the Bargmann index M/2. By 



making use of the disentangling theorem of su(l, 1) algebra [[H]] we arrive at the displacement 
operator formalism of NBS 



\r]e 10 ; M)~ = exp [C c /C+ - [| 0> , ( c = e jy arctanhr/. (3.5) 

Eq.( p.5[ ) is nothing but the Perelomov's coherent state of su(l, 1) via its HP realisation. 
The Perelomov's coherent states admit the ladder operator form [IS|. To see this, we 



differentiate both (|3.3|) and ( |3.5|) with respect to \(\ and equate the results. We have 



e~ l& /C_ - v e iy JC+\ \r]e M ; M)~ = M^e 10 ; M)~. (3.6) 

This ladder operator form is obviously compatible with the limit results. In fact, in the 
limits r] — > and rj — > 0, M —>■ oo with Mrf = a 2 , Eq.( [3.6|) reduces to 



6|0;M)~ = 0, 6|0;oo) _ = ae l<, |0;oo) _ , (3.7) 
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which mean |0; M)~ and |0; oo)~ are the vacuum and coherent states respectively. 

Finally we point out that the NBS can also be regarded as the density dependent anni- 
hilation operator coherent states, namely, 

E M \^ e - M)~ = e id rj\rje i6 ; M>", E M = b-^=J== (3.8) 

when M > 2. The operator E± is not well-defined in the vacuum state. If we require 

J5r||0> = 0, (3.9) 

then it is obvious that E± = E~ . This connection of SG phase operator with Perelomov's 
coherent states with Bargmann index 1/2 was well-known. 

4 Nonclassical effects 

4.1 Photon statistics 

Let us first examine if the NBS is of sub-Poissonian statistics. Using Eq. fl2.5|) , it is easy to 
calculate the averages (N), (N 2 ) and the fluctuation (AN 2 ) of the photon number 

Then Mandel's Q-factor characterising sub-Poissonian (if Q < 0) distribution is obtained as 

W (AT) 1 - 1) 2 1 ' 

So the field in NBS is super-Poissonian, not sub-Poissonian, except for the (vacuum and 
coherent state) limit rj — ► (Q — > 0, Poissonian statistics). 

Since the occurrence of antibunching and sub-Poissonian are concomitant for single and 
time independent field, as in the present case, the field in NBS is bunching, not antibunching. 
In fact, from the following second-order correlation function g 2 (0) 

2/ . (6+6+66) 1 , x 

9 (0 > = W = 1 + m > (4 - 3 > 

we can arrive at the same conclusion. 

4.2 Squeezing effect 

Let us evaluate the variances (Ax 2 ) = (x 2 ) — (x) 2 and (Ap 2 ) = (p 2 ) — (p) 2 of the quadrature 
operators x (coordinate) and p (momentum) defined by 

x = ^ + b), p=-L(&t_ 6 ). ( 4. 4) 



From the following relation 



b«\r]e i6 ;M) 



r/e 



to 



M{M + !)••• (M + k — 1) \r,e ie - M + k)~ 



(4.5) 



we have 
(Ap 2 ) - 



1 , 2 (1-^ M » 2 „(M + n-l)! 



2 (M-l)! 
-2sin 2 6»r] 2 (l - r? 



n=0 



2\2M 



III 



E 

n=0 



M + n-1 
n 



y/M + n (cos(26)VM + n + 1 - VM + n) 

(4.6) 



7? 2r VM + n 



(Ax 2 ) 



1 rfM + cos(2% 2 (l -r/ 2 ) M x 



(M-l)! 



y2 v 2n(M + n-l)\ 



-2 cos 2 9 rf {I - rj 



2\2M 



E 

.n=0 



n=0 

M + n-1 



n! 



(M + n + l)(M + n) 



r/ 2n VM + n 



(4.7) 



In the derivation of Eq.( [4.6| ), we have used the identity 
Mr] 2 rf(l-7] 2 ) M Ml 



(M-l)! (l-ri 2 ) 



M+l 



V 2 (l~V 2 ) M f^^M + n) 



(M-l)! 



n=0 



nl 



r) 2 (l-rj 2 ) M « 2n (M + n-1)! 



(M-l) 



E^ 

n=0 



n! 



VM + nVM + 



n. 



(4.8) 



Let us analyse (Ap 2 ) in more detail. First consider the case 8 = (the initial time). In 
this case Eq. ( |4.6| ) is simplified as 



?7 2 (1 — ?7 2 ) 



2\M oo 



g ^2n (M + nJ)! - / v/M + n + 1 _ - 



n) . (4.9) 



- (M-l)! n=0 

Since every term in the infinite series in Eq. (|4.9|) is positive, the infinite sum is positive. So 
we always have (Ap 2 ) < 1/2. This means that the quadrature p is squeezed. Fig.l shows 
how P = (Ap 2 ) depends on 77 and M. We have chosen < rf < 0.99. Those plots show 
that: 

(1) . Dependence on rj 2 . It is found that the variance (Ap 2 ) is an increasing function 
of rj 2 , namely, the larger rj 2 , the stronger the squeezing effect. Note that (Ap 2 ) > since 
(Ap 2 ) = (p 2 ) and p 2 is a positive definite hermitian operator. 

(2) . Dependence on M. We have chosen M=l, 5, 50. The squeezing of p for larger M 
is stronger than that for small M. When rj 2 is small (close to 0) or large (close to 1), the 
difference is very small. While when 11 2 is around 1/2, the difference is larger. However, 
squeezing is not so sensitive to M. The case M = 5 and M = 50 are almost same as showed 
in Fig.l. 

Furthermore, due to the uncertainty relation (Ax 2 )(Ap 2 ) > 1/4, the variance (Ax 2 ) > 
1/2 when 9 = 0. So the quadrature x is not squeezed. 
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When 9 ^ 0, it is easy to see that (Ap 2 ) is a 7r-periodic function with respect to 9 and 
symmetric with respect to 9 = ir/2. To investigate the effect of 9 to (Ap 2 ), we plot the 
(Ap 2 ) as a function of 9 for different rj 2 values in Fig. 2 (we choose M = 1 for simplicity). 
We find that (1) (Ap 2 ) becomes larger when 9 creasing. It first reaches 1/2 and then reaches 
the maximal value when 9 = tt/2. Then it symmetrically decreases until 9 = tt. In some 
region around 9 = n/2 (dependent on r/ 2 ), the squeezing disappears. (2) Small rf violates 
the squeezing slightly while large rj 2 destroys the squeezing strongly. 

5 Generation of NBS 

The displacement operator formalism suggests that the NBS can be generated by the non- 



degenerate parameter amplifier described by the Hamiltonian |]T9 

H = H + xi(a[ale~ 2luJt - a x a 2 e lwt ), H = u^ai + u 2 a\a 2 , (5.1) 

where a± and uo\ (a 2 and u 2 ) are the annihilation operator and frequency for the signal 
(idler) mode. Frequencies iO\ and oj 2 sum to the pump frequency, 2ui = uj\ + uj 2 . The 
coupling constant x is proportional to the second-order susceptibility of the medium and to 
the amplitude of the pump. The unitary time evolution operator in the interaction picture 
is 

U{t) = e iHot e x * (a i a 2~ aia2 ) e ~ iH o t _ (5-2) 

Suppose that the system is initially prepared in the state |0,M) = ||0). Since the photons 
are created or annihilated in pairs, we can restrict ourselves in the subspace 

\\n) = \n,n + M), n = 0, 1,2, (5.3) 

which is isomorphic to the single-mode Fock space. Then at any time t the system is in 



tf(*)l|o> = E 

n=0 L 



U + M — 1 \ , ^ _ 2 _ ,\Af/j.„_i. . ,^2n 



77 



2 



e l2un \\n). (5.4) 



l-tanh 2 xt) (tanhxt) 
Identifying rj = tanh^t and 9 = 2ut, we obtain the NBS. 

6 Conclusion 

In this letter we have introduced the negative binomial states and studied their nonclassical 
properties. The following table shows some properties of NBS with the special emphasis on 
the comparison with the binomial states: 
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NBS 


BS 


Limit States 


Coherent and SG phase states 


Coherent and number states 


Related algebra 


su(l,l) 


su(2) 


Essence 


Perelomov's coherent states 


Atomic coherent states 


Sub-Poissonian 


No 


Unconditional 


Antibunching 


No 


Unconditional 


Squeezing effect 
(0 = 0) 


Unconditional for p, 
no squeezing for x 


Limited ranges of parameters 
for x, no squeezing for p 



We see that the nonclassical properties of the NBS and BS are complementary. 

As further work we shall generalize the notations of NBS in this letter to the negative 
multinomial states with negative multinomial distribution as their photon distribution. This 
generalisation concerns the multi-mode radiation field and should exhibit some more fruitful 
nonclassical properties like correlation between different modes. It is also a good challenges 
to study the interaction of radiation field in the NBS with the atoms. 
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g -]_. Gray line: M=l 
Black line : M=5 
Dashed line: M=50 

° 0.2 0.4 0.6 0.8 1 

T|**2 

Figure 1: Variance P = (Ap 2 ) as a function of rf (= rj * *2) for 6 = (initial time) and 
M=l (gray line), 5 (black line) and 50 (dashed line). 




e 

Figure 2: Variance P = (Ap 2 ) as a function of 6 for different r^ 2 , 0.2 (a), 0.5 (b), 0.7 (c) 
and 0.9 (d). In all the cases M — 1. The dashed line corresponds to rj = or P = 0.5. 
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